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Abstract 

Oscillons are long-lived, localized excitations of nonlinear scalar fields which may be copiously 
produced during preheating after inflation, leading to a possible oscillon-dominated phase in the 
early Universe. We investigate the stochastic gravitational wave background associated with an 
oscillon-dominated phase. An isolated oscillon is spherically symmetric and does not radiate grav- 
itational waves, and we show that the flux of gravitational radiation generated between oscillons 
is also small. However, a significant stochastic gravitational wave background may be generated 
during preheating itself (i.e, when oscillons are forming), and in this case the characteristic size of 
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5^ the oscillons is imprinted on the gravitational wave power spectrum, which has multiple, distinct 

peaks. 
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I. INTRODUCTION 

After inflation ends, the energy of the coherent homogeneous infiaton field must be re- 
leased to reheat the Universe, giving way to the hot Big Bang. This may occur via pre- 
heating [IHZ1 where particles are explosively produced with non-zero momenta, followed by 
thermalization. In some cases, high-momentum quanta of the infiaton field can be copi- 
ously produced before other particles [8hTT]. The particle production is seeded by initially 
minute quantum field fluctuations which are exponentially enhanced via parametric reso- 
nance, driven by the homogeneous, oscillating infiaton field. This process is nonlinear and 
nonperturbative. Fortunately, as the particle number increases exponentially the occupa- 
tion numbers of higher momentum states grow rapidly, and the quantum fields can then be 
well-described by stochastic classical fields [SJ. 

Many interesting non-perturbative field theory phenomena may take place during pre- 
heating, including formation of topological or non-topological solitons p21 [13] as well as 
quasi-solitons such as oscillons [HI IHJ [HH3I]. Oscillons are spatially localized, temporally 
oscillating configurations of a nonlinear scalar field theory. First described by Bogolubsky 
and Makhankov |15j . advances in numerical methods have allowed these fascinating objects 
to be studied in great detail [T6T431] . They are not strictly stable, even in the classical limit, 
as they are not associated with conserved charges (topological or non-topological), but are 
long-lived due to a trade off in the oscillon potential between the dissipative quadratic terms 
and an attractive force generated by the nonlinear terms. 

A number of models have been described in which massive oscillons are generated during 
preheating [8HTT 1 H H 135] . In particular, Ref. [8] showed that in a large class of well-motivated 
inflationary scenarios oscillons may account for the majority of the energy density in the 
post-inflationary Universe. In these scenarios, the Universe undergoes a transient (effective) 
matter-dominated phase. For a given inflationary model, a transient matter-dominated 
phase modifies the effective spectral index of the inflationary perturbation spectrum relative 
to that produced via immediate thermalization [361 EZ] > although this effect degenerates with 
other mechanisms that cause the expansion history of the early Universe to depart from a 
strict radiation-dominated scenario. 



In this paper, we examine possible stochastic gravitational wave signals generated during 
this oscillon preheating scenario, to see whether there are any characteristic features that 
are - even in principle - observable today. Gravitational waves, once generated, propagate 
freely even in the very hot, dense early Universe, and are already known to be generated 
during resonance [3BHJ7]. Typically, gravitational waves generated by preheating after "high 
scale" inflation are found in the MHz to GHz range, well beyond the scope of current or 
planned interferometric detectors. However, their importance as a potential probe of the 
post-inflationary Universe and inflation models has motivated the study of the possible 
stochastic backgrounds generated during preheating and provided impetus for the develop- 
ment of novel detector technologies. 

In what follows, we will show that the characteristic scale of the oscillons is imprinted on 
the gravitational wave background associated with the resonant phase that generates them, 
giving rise to power spectra with organized broad peaks. Conversely, once initial transients 
decay, oscillons are rather spherically symmetric, so an isolated oscillon does not radiate 
gravitational waves. We also compute the gravitational wave emission generated between 
multiple oscillons, and show that this is very small, so an oscillon-dominated phase itself 
does not generate significant gravitational waves. However, the amount of gravitational 
wave energy produced when the oscillons are forming is typical of that of a usual preheating 
scenario. 

The outline of the paper is as follows. In Section |TTJ we briefly review oscillon dynamics 



and introduce the oscillon preheating scenario. In Section III we analytically estimate the 



gravitational wave production from the oscillon-dominated phase and show that the pro- 



duction is highly suppressed when oscillons are properly formed. In Section |IV| we describe 
our numerical algorithm and present the results of our lattice simulations, and we conclude 
in Section IVl 



II. OSCILLONS AND OSCILLON PREHEATING 

Oscillons can form in a variety of models (see e.g. [8j |20j |31j [35]) , but given that an 
oscillon-dominated Universe can follow inflation driven by a single minimally-coupled scalar 
field, we restrict our attention to this scenario [SJ. Oscillons can form when the scalar 
potential "opens up" [31]. Suppose the potential has the form 
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n0) = T 2 + Ki(0), (l) 

with a minimum at V(<fi = 0) = l and is symmetric under — > — 0, then the 'opening up' 
condition simply means [Mj 

V n \((f>) < for some region of 0. (2) 

To understand this condition, we will need an approximate analytic description of oscillon 
physics. For oscillons much smaller than the characteristic spacetime curvature, the Klein- 
Gordon equation is simply 

4> - V 2 + m 2 + ^(0) = 0. (3) 



1 A minimum at non-zero 4> may be shifted to the origin by a field redefinition. 



An oscillon usually has a dominant oscillating frequency which we write as 

4>(x) = $(x) cos(wt), (4) 

where $(x) is a positive localized function that falls off quickly far from the center. With this 
ansatz, far from the center of the oscillon the field equation reduces to (m 2 — w 2 )$ ~ V 2 $. 
Since V 2 $ has to be small but positive in this region, the angular frequency satisfies 

oj ~ m and u < m, (5) 

showing that the mass parameter sets the dominant oscillation frequency. In the central 
region, the nonlinear term V^ is important and V 2 $ is negative. In combination with Eq. (15]) 
this leads to — oj 2 $ — V 2 $ + m 2 $ > 0. Since Eq. 0) should be satisfied for any t, we get that 
V^[(<I>) = — (— u 2 Q— V 2 $+m 2 $) is negative in the oscillon center if cos(mt) = 1, and V^(— $) 
is positive it cos(mt) = — 1. Integrating over the field from to 0, we then obtain condition 
Q2J) . Physically, if the interaction term is shallower than the mass term (e.g. —(f) 4 ), it induces 
an attractive force between "particles", which can balance the dissipative quadratic terms, 
and rendering oscillon formation energetically preferred. 

Our ansatz Q captures the dominant oscillating mode of the oscillon. However, an 
oscillon also has higher order harmonics; e.g. cos(2mt), cos(3m£), ■ • • [281 [33]. We can see 
this nicely in the small amplitude limit [28J where the smooth configuration scale of the 
oscillon is much larger than its natural oscillation period, i.e., 

d d , 

— ~ m ~ — — - and <p ~ em, (6) 

at eax % 

e being a small parameter. To count orders, we expand the field as 



[X) = 

n=l 



T,<?Mx), (7) 



and redefine the spatial coordinates x % = ex 1 . Assuming V n \((j)) = J2 a >2 9a4 )a / a (9a being 
coupling constants), the equations of motion for the modes n are 

(f) n + m 2 (j) n - V?0 n _2 + ^J ^J ga<f>m<f>n 2 -~<f>n a -i = ° (8) 

a>2 ni+n,2...+n a -i=n 

where n goes from 1 to oo and we define 0_i = <po = 0. We can solve this set of equations 
perturbatively, and the first two equations are 

0i + m 2 0! = 0, (9) 

2 + m 2 2 + g^l = 0. (10) 

Eq. ^ gives 0i = $i(x) cos(mt), recovering our ansatz Q. (To be fully rigorous here, 
we should have allowed for a spatially dependent phase in the cosine, but the boundedness 
of the whole solution forces it to be constant [28], and we have set it to zero.) Inserting 



our solution for 1; Eq. (10) reduces to a forced harmonic oscillator, with a forcing angular 



frequency 2m, with the solution 



$21 (x) cos(mt) + $22 (x) sin(mt) + 93 ^/ [cos(2mt) - 3]. (11) 

6m 2 



By induction we see that the higher order equations are all forced oscillators, and at each 
step in this hierarchy the (largest) forcing angular frequency increases by m, which in turn 
increases the (largest) angular frequency for <f> n by m. If the potential is symmetric under 
— > —(f) around the minimum, the odd order coupling constants vanish ((fea+i — 0). From 



Eq. ( 11 ), we see that then the second harmonic cos(2m£) vanishes and this can be generalized 
to higher orders, so higher harmonics with the angular frequency 2nm all vanish. The 
existence and properties of these various harmonics away from the small amplitude limit 
have been investigated numerically in [33J , with results consistent with this small amplitude 
analysis. 

From an inflationary perspective, the "opening up" condition (|2| is relatively easy to 
satisfy as current data disfavor the quadratic potential relative to flatter potentials [37, 48- 
[30] . For a consistent description of preheating with a single scalar, we must also assume 
the inflaton's self-couplings, at least initially, are more important than its couplings to other 
fields. Inspired by a number of string and supergravity inflation models [5TT456] , a family of 
potentials have been proposed that yields an oscillon preheating scenario whilst satisfying 
current observational constraints [8] 

m 2 M 2 



V&) 



la 



1 + W 



(12) 



where < a < 0.9 and M < 0.05M P (M 2 P = l/8nG). In particular, a = 1/2 is the axion 
monodromy model [511 E2]- When <p/M is large, the potential can be approximated by 
V(4>) ~ m 2 M 2 ~ 2a (f) 2a /2a. On the other hand, when <j)/M is small, we have V(4>) ~ m 2 (j) 2 /2. 
For this potential very massive oscillons can be copiously produced after inflation, giving rise 
to an oscillon-dominated epoch before thermalization. In this paper we ask whether oscillon 
preheating leaves any "footprints" in today's Universe, particularly in the gravitational wave 
background it generates. 



III. GRAVITATIONAL WAVES: ANALYTICAL STUDY 

There are three possible ways where gravitational waves can arise as a consequence of an 
oscillon-dominated phase in the early Universe - during their formation, during the oscillon 
dominated phase itself, and during their decay and the subsequent thermalization of the 
Universe. We cannot address the last issue, since this process is beyond the scope of the 
model we are analyzing. The first process will be driven by the strong nonlinear dynamics 
of parametric resonance and will be investigated numerically in the following section. In 
this Section we give a semi-analytical treatment of gravitational wave production in the 
oscillon-dominated phase itself. 

In simulations of oscillon formation, the resulting oscillons are essentially at rest with 
respect to the background spacetime. This makes sense: we can view an oscillon as a 
collection of a scalar particles (i.e. , a system with a large number of degrees of freedom) which 
are generated from a homogeneous field configuration with a vanishing total momentum, so 
the overall momentum of the generated oscillon should be highly suppressed. Consequently, 
as a good approximate description, we can assume that oscillons are not moving. Secondly, 
oscillons are formed by causal interactions, well within the Hubble horizon, so we estimate 
the production of gravitational waves from oscillons in Minkowski space, where we can 
use standard gravitational wave formulae, summarized in Appendix |Aj Thirdly, gradient 
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FIG. 1: A late time snap shot of the energy density in oscillon preheating. The model is that 



of Eq. (12) with a = 1/2 and M = O.OlMp. The box size is L = 50/m and the energy density 



isosurface is taken at a value 5 times the average energy density. 



energy is minimized by spherically symmetric profiles, and numerical simulations show that 
the oscillons approach a symmetric state as they evolve (see Figure (0])). Consequently we 
assume spherically symmetric oscillons. (In Appendix pi we show that our results also apply 
for the case of small deviations from perfect spherical symmetry.) Finally, the nonlinear 
potential makes it difficult to analytically describe the spatial profile of the oscillon and we 
will assume a Gaussian spatial profile for convenience. 

Consider a Minkowski space homogeneously distributed with the idealized oscillons de- 
scribed above. According to Birkhoff's theorem, there is no gravitational wave emission 
from a spherically symmetric source, oscillating or not. As a consistency check, one can 
substitute the ansatz </>(£, x) oc e~ m ° x I 2 cos(/it) (where m and \i characterize respectively 



the width and angular frequency of one harmonic of the oscillon) into Eq. (A5) and find 
that the emitted power in gravitational waves indeed vanishes. 

Gravitational wave production can still be driven by pair-wise and higher order "interac- 



tions" between oscillons. Consider the following ansatz 



4>(t, 



Y^Aae-^*-^ 2 cos(/i£ + ip a ) 



(13) 



where A a is the amplitude of the center of the oscillon, d a is the spatial position of the oscillon 
and ip a the initial phase. The resulting gravitational waves have angular frequency 2/i, since 

;0 ~ COS 2 (/it + if a), 



schematically the perturbative Einstein equation is Ohij ~ T^ ~ d. 

where hy is the perturbative metric. Substituting this ansatz into the formula (A5), after 

some tedious Fourier transforms and algebra we get the emitted power 

dP 2M _ dP 2 dP 4 



with 



dfl 

dP 4 
dft 



dn dn dn 



2 s~i 2 

" T Y,( A aA)\fid ab )\m d ab ) 2 sm\±,d ab )e^ mldlb ~^ 
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a>b 



n 2 G 
l6~ 



^2 A a A b A c A d (fxd ab ) 2 (m d cd ) 2 S abcd (±, d, ip)e 



1 2 fd 2 , d 2 \ » 



(a>6)^(c>d) 



where 5^ a>b implies summation over a and 6 that satisfy a > b and J^/ a 



(a>6)^(c>d) 

6, c and d that satisfy a > b, c > d and (a, 6) 7^ (c, d). We have introduced x 
d ab = d a - d b} ip ac = ip a - ¥c, and 



(14) 

(15) 
(16) 

is over a, 

x/|x|, 



Sabcd (X, d, If) 



[cos(d afe , d cd ) - cos(x, d ab ) cos(x, d cd )] 2 - - sin 2 (x, d ab ) sin 2 (x, d cd ) 

■ C0S(— =X ■ (jidac + Hd bd ) + ip ac + lf bd ). 



(17) 



where sin(x, d a b), denotes the sine of the angle between x and d ab . dP2/dQ is the contribution 
from pair-wise oscillons and dP^/dQ comes from the four-oscillon contribution, and because 
of the structure T^ ~ di<frdj<f>, there is no other contributions. 

The dominant frequency of oscillons is slightly smaller than m (the mass parameter in 



the potential) and most of the energy-density is localized within a radius r 



m 



of the 



oscillon core, so as a good estimate we can set m = \x = m in the rest of this Section. (The 
results are similar if we consider the higher harmonics \x = 2m, 3m, ....) 

We can see that either the two-oscillon or the four-oscillon contribution is exponentially 
suppressed (i.e. m 2 d 2 b ^> 1), since in the oscillon preheating scenario the separation of 



oscillons are usually 



> 



m 



For the two-oscillon contribution, this simply means that the 



leading contribution is from immediately neighboring oscillon pairs. Thus the power density 



of gravitational wave emission (the power of emission (14) divided by the space volume) 
from dP 2 /dfl is highly suppressed. 

The suppression of the four-oscillon contribution is less straightforward to see, even if 
we intuitively expect it to be small. The exponential suppression above is by (md ab ) 2 
and (md cd ) 2 . While these suppressions restrict attention to small {md ab ) 2 and (md cd ) 2 
contributions, they do not suppress large separations between oscillon a and oscillon c and 
between oscillon b and oscillon d. That is, in the summation ^( a>b )-^( c>rf ), it might appear 



that we have to consider a large number of four-oscillon contributions (the ones with small 
(rndab) 2 and (md C( j) 2 but (arbitrary) large (md ac ) 2 and (md^) 2 )- If this were the case, 
the four-oscillon contribution to the power density of gravitational wave emission would 
diverge, contradicting the usual principle that energy is extensive and can be considered as 
a sum of those of the sub-volumes. To see how locality is respected, focus on the factor 
cos(x • (md ac + mdbd)/V2 + (fac + fbd), which is the only place where d ac and du enter the 
formula. When (md a (,) 2 and {md C( i) 2 are small and d ac and d^ are large, d ac and d bd have 
to be almost parallel, in which case 

cos(— =x ■ (md ac + md hd ) + ip ac + ip bd ) ^ cos(V^|md ac | cos 9 + <p ac + <p bd ), (18) 

9 being the angle between x and d ac . When |md ac | is large this cosine factor oscillates 
rapidly about zero when cos# varies (Note that dQ = dcos#d0.). Physically, the energy 
emission is obtained by integrating a finite solid angle. For a given solid angle, large enough 
|md ac | this cosine oscillation suppresses any net contribution. Thus only configurations with 
four adjacent oscillons can contribute and these, as mentioned above, are highly suppressed. 
This concludes that gravitational wave production is negligible during the oscillon dominated 
phase itself. 



IV. GRAVITATIONAL WAVES: NUMERICAL SIMULATIONS 

In this Section, we compute the gravitational wave production during the resonance and 
oscillon-formation phase. Our numerical code is based on PSpectRe [57] , a pseudo-spectral 
code for simulating the evolution of coupled scalar fields on an expanding background, 
upgraded with the capacity of computing gravitational wave power spectra. The initial in- 
homogeneous seeds of the scalar field are generated similar to DEFROST |58j. The scalar 
field evolution is computed using a second-order-in-time symplectic integration scheme with 
a fixed time step. The transverse-traceless part of the first order perturbative metric is 
computed using the method described in [42J - with the Fourier modes of the metric per- 
turbations (and their sources in T^) evolved in Fourier space, via a fourth-order-in-time 
Runge-Kutta scheme. To first order in the metric perturbations, the entire combined sys- 
tem is self-consistent: the background evolution is self-consistently computed from the scalar 
field and the evolution equations used to evolve the metric perturbations are valid at all scales 
and do not assume any particular expansion history. 



The results presented here are for the model of Eq. (12) with a = 1/2 (i.e., the axion 
monodromy model) and M = O.OlMp. The oscillon mass m is determined by matching the 
amplitude of the temperature anisotropies in the cosmic microwave background 2 



1 I m \ I M 



A 2 =^— — — (AaAf) 1+a = 2.4 x 1(T 9 (19) 

R 96n 2 a 3 \M P ) \M P ) K ' K ' 

where J\f, the number of e-folds before the end of inflation, is set to 55 here. We have 
checked the stability of our results against changes in lattice and box sizes as well as physical 



2 We used the amplitude of the power spectrum obtained from the WMAP 7-year analysis [48 , but our 
results are not strongly dependent on this parameter. 
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FIG. 2: Evolution of the gravitational wave power spectrum. The lattice is 256 3 , the box size is 
L = 50/m and the evolution time shown is from t = 80/m to t = 240/m. There are (at least) 
four distinct stages: (top-left) - steady "pumping" in low momentum modes, linear perturbation; 
(top-right) - rapid growth, oscillons forming; (bottom-left) - oscillons stabilizing, peaks and troughs 
becoming significant; (bottom-right) - long stable stage, very slow growth. 



parameters such as M. The results presented here are derived from simulations on a 256 3 
lattice with box size L = 50/m (the initial Hubble horizon being about 270/m), unless 
otherwise stated. 

Figure pi) shows the time evolution of the gravitational wave power spectrum, plotted 
using the present-day frequency and fractional energy density in gravitational waves, as- 
suming instantaneous thermalization at the instant the spectrum is computed (see e.g. |Hj 
for the conversion.) If there is a long oscillon-dominated phase, this last assumption will 
need to be revisited, as we discuss below. The power spectrum curve grows upwards with 
time, which is labelled in units of 1/m. In what follows, we describe in succession four 
stages that can be easily identified. A significant production of gravitational waves begins 
at around t = 80/m, with an initial steady growth, characterized by the "pumping" of the 
low frequency modes (top- left graph). Here fluctuations of the inflaton field are small com- 
pared to the zero-mode oscillation. During this stage the total power fluctuates as newly 
created gravitational waves can constructively and destructively interfere with those that 
are already there, although the overall trend is toward greater power. The second stage is 
rapid growth (top-right), as oscillons emerge and the zero-mode oscillation of the inflaton 
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FIG. 3: The energy densities of the same slice of the box at two different times. The model is that 
of Eq. |l2} with a = 1/2 and M = 0.01M P . The lattice is 256 3 and the box size is L = 50/m. 



The left is taken at the late time of the rapid growth stage and the right is taken at the early time 
of the stabilizing stage. We can see that at the onset of the stabilizing stage oscillons become very 
massive, dominating the energy density of the Universe. 



decays rapidly. Interesting peaks and troughs starts to emerge in the power spectrum at the 
later time of this stage. The rapid growth stage is followed by a stabilization, as oscillons 
further drain energy from the environment, with existing oscillons becoming more massive 
and new smaller ones emerging (bottom-left graph). Figure (J3J) shows the energy density 
as a function of position on a spatial slice at two close times, one taken at the late time 
of the rapid growth stage and the other taken at the early time of the stabilizing stage. 
We can see that oscillons in the stabilizing stage are much more massive than in the rapid 
growth stage. More interestingly, back to Figure ([2]), the structure of peaks and troughs 
becomes more evident at the stabilizing stage. Finally, we have a long stable stage, when 
oscillon formation is complete and the gravitational wave power spectrum grows very slowly 
(bottom- right graph). 

In order to better investigate the peaks in the spectrum we repeat these calculations in 
a smaller box (25/m) to resolve higher frequencies, with results shown in Figure (El). As we 
discussed in Section |TTJ the oscillon has a dominant angular frequency of u ~ m, as well as 
higher order harmonics with angular frequencies 3m, 5m, 7m, .... The even order harmonics 
are absent, because the model given by Eq. (12) is symmetric under — > —<f). Interestingly, 



we see the same structure in the gravitational wave power spectrum - the spacing of the 
troughs has the same pattern as the oscillon harmonics, as also shown in Figure (El). 

We have shown in Section III that the gravitational wave emission from spherically sym- 
metric oscillons is highly suppressed. Therefore, the fact that these troughs and peaks be- 
gin to appear when oscillons are becoming massive and increasingly spherically symmetric 
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FIG. 4: Peaks and troughs of the gravitational wave power spectrum. The model is that of Eq. (12) 
with a = 1/2 (the axion monodromy model) and M = O.OlMp. The lattice is 256 3 , plotted at 
t = 240/m. The box size is L = 25/m, instead of L = 50/m, so as to resolve more high momentum 
modes. The vertical lines correspond to the gravitational wave frequencies associated with the 
different harmonics of the oscillon, which are twice the frequencies of the oscillon harmonics. Even 
order oscillon harmonics (2m/ ir, 4m/ tt, ..., in terms of the gravitational wave frequency; the dotted 
lines) are absent in our model, because of the symmetry <fi — > —<fi. The matches of the vertical 
lines with the troughs indicate that oscillons are suppressing the gravitational wave production at 
those frequencies. 



seems to suggest a simple explanation: the troughs are forming because spherically symmet- 
ric oscillons suppress the power growth at those bands. As we have mentioned in Section 



IITJ the angular frequencies of the gravitational waves these different harmonics excite are 
2m, 6m, 10m, 14m, ... respectively. In Figure (El), we have plotted the gravitational wave 
frequencies (the angular frequencies divided by 2ir) associated with the different harmonics 
of the oscillon - they match the troughs rather well. 

Oscillons are quasi-stable but our gravitational wave spectra are plotted using a trans- 
fer function which assumes instantaneous thermalization. During the matter-like oscillon- 
dominated phase the gravitational waves will redshift relative to the overall density of the 
Universe, reducing their energy density, but their present-day frequency is also reduced, since 
the comoving horizon size is growing more slowly. Since the overall duration of the oscillon 
dominated phase is unknown (it depends on couplings between the irmaton and other fields) 
we can only compute this correction parametrically. If the Universe grows M OSC i e-folds 
larger during the oscillon-dominated phase, gravitational wave frequencies are redshifted by 
a factor e A/ ° sc */ 4 while the power spectrum amplitude is diluted by a factor of e osci . Also, a 
modified expansion history affects the evaluation of M which determines the value of m via 
Eq. (19) and thus affects the gravitational wave power spectrum. For example, if N sa ~ 7 



and assume the gravitational wave production is not significant in the decay of oscillons, 



11 



considering all these effects, we would have f^ tcal ~ 10 8 Hz and h 2 VL t ^ %cal ~ 10 14 . 

V. CONCLUSIONS 

We have studied the stochastic gravitational wave background associated with an oscillon- 
dominated phase in the early Universe. We see significant gravitational wave production 
as oscillons are being formed, and the resulting spectrum has a specific and characteristic 
structure, with multiple peaks and troughs as shown in Figure Q. Further, the frequencies 
associated with the peaks and troughs are matches to those that contribute to the internal 
dynamics of the oscillon. Conversely, when oscillons are properly formed and dominating 
the Universe, the gravitational wave production is not significant. Both the appearance of 
multiple troughs in the power spectrum and the lack of gravitational energy emission in the 
oscillon-dominated phase itself seems to stem from the fact that oscillons, after formation, 
morphologically evolve towards spherical symmetry. 

Generically, considering stochastic gravitational backgrounds produced at the end of in- 
flation leads us into a somewhat paradoxical position. Intellectually, these are among the 
most fascinating topics in cosmology, and provide a potential window into the Universe at 
the end of inflation. From a practical perspective, however, they occur at frequencies far 
beyond the capabilities of present-day gravitational wave detectors, or even any reasonable 
extrapolation of interferometer-based technologies. Consequently, detecting a stochastic 
gravitational background generated at the end of "high scale" inflation will require some 
new technology, and the signal discussed here is no exception. 

This analysis deepens our understanding of the phenomenology of monodromy inflation 
and its generalizations. From a theoretical perspective, monodromy inflation is one of the 
most solidly-motivated inflationary scenarios [5T| 152"] and a good match with current data 
[55] . In addition to a possible oscillon-dominated phase, monodromy inflation is also as- 
sociated with a possible modulated power spectrum, a possibility which can be usefully 
constrained with modern astrophysical datasets [60J and would be correlated with a sub- 
stantial 3-point function [SI) . Consequently, it is worthwhile to explore the properties of 
this model. The rich phenomenology of monodromy inflation also underlines the impor- 
tance of carefully considering the full consequences of inflationary models, in order to better 
understand - and eventually test - their predictions. 

The appearance of multiple troughs, as we argued, is due to formation of spherically 
symmetric oscillons and thus presumably is quite universal in models where oscillons can 
form and dominate the Universe. It would be interesting to see whether oscillon preheating 
can arise in any "low scale" inflation models and thus leads to signals detectable in the near 
future. 
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Appendix A: Formulae of the gravitational wave power spectrum 

Given a source with energy momentum tensor T^ u (t, x), the gravitational energy emitted 
per solid angle is given by [62J 

/■oo 

AiiMk) / co 2 duT l \u,k)(T lm (u,k)y. (Al) 

Jo 



dE G 

dQ ~ 2^' 



where k = uk = ux = wx/|x|, the transverse-traceless projector is defined as 

l\-ij,lin\K) = OuOjm 2KjrC m du + — fcjfcjfc;fc m ~OijOi m + —OijKiK m + — 0/ m fcjfcj (A2J 

and the Fourier transform of the energy momentum tensor is 

T -(w, k) = / d 4 xe^- lk - x T 4J (t, x). (A3) 



Since Aij s i m (k) projects out the trace part of 7#(u;, k), we can drop the manifest trace term 
5ijC in Tij(t, x) to simplify the calculation. Note that this formula has taken into account the 
conservation of the energy-momentum tensor and thus the scalar's equation of motion, as 
for a scalar its equation of motion is equivalent to the conservation of its energy-momentum 
tensor. 

If there are only some isolated frequencies un{N = 1,2,...) in the energy-momentum 
tensor 

T^ u (u, k) = ^ T^iu, k)27T(5(u; - u N ), (A4) 

JV 



then we can rewrite Eq. (Al) as 



dP y, GtU 2 

N 



^ == i^ —^k ihlm (k N )T\uj N , k N )(f lm (u N , k N ))* (A5) 



where again kjy = watx and P = E/time is the power of the energy emission. Because 
of the projection Ayy m (A;), once the field profile is known, the potential of the scalar field, 
which only appears as a trace in T^-, is not needed in the estimate of the gravitational wave 
production. The potential information in a sense is already encoded in the field profile. 

Appendix B: Deviation from spherical symmetry 

Oscillons are most likely not perfectly spherical, especially at formation. Here we show 
that the suppression is still significant if the deviation from the spherical symmetry is not 
dramatic. 

We can imagine an infinite number of oscillons homogeneously distributed across the 
space, with the two-oscillon and four-oscillon contributions being exponentially suppressed, 
as with the spherically symmetric case. So, here we focus on the contribution from a single 
oscillon and adopt the following ansatz 

4>(t, x) = Ae-^ m " x2+m y y2+m " z2) cos(/it), (Bl) 
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FIG. 5: The function F(a)/a 2 . m x = m y = fJ,,m z = afi, where a encodes the extent of the ellipse 
of the oscillating oscillon. For a = 1, we have -F(l) = 0, which means the oscillon is spherically 
symmetric and around which the gravitational energy emission is suppressed. 

where m x , m y and m z are in general different, thereby breaking spherical symmetry. Again, 
after some Fourier transforms and algebra, we obtain the power per solid angle as 



in /i4,,2 



dP 2fM 2tt 2 GA*ii 



dfi (m x m y m z ) 



, J~ \J>vii •£{ 



(B2) 



where 



£ 



mf(l- 2xi 



+ £ 



2 »2 



£ 



m 2 (x- - 1] 



VEi 



The integration J dfiJ-^m^Xj) is difficult to do analytically for general m^, m y , m z and /i. 
As a representative example, we can consider the case that the oscillon is elliptical, that is, 

m x = m y = fi, m z = aft, (B4) 

where the dimensionless a encodes the extent of the ellipse. For this case, we have 



yw^H^w 



(B5) 



with 



7T 



F(a) = ^e-^V(-Ila 2 + 



^3/2 erf (^/4/cr 2 -4 

^(83a 4 -144a 2 + 64) V 



>/4/cr 2 -4 



(B6) 



where erf () is the error function. The integrated power of the gravitational wave emission 
for an elliptical oscillon is given 



2p 



in /i4 



2tt 2 GA 



F(a) 



a" 



(B7) 



where F(a)/a 2 is numerically given by Fig. rt5J). We can see that the curve is very flat 
around a — 1, the spherically symmetric case. So suppression in the gravitational wave 
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production is still significant even if the oscillon shape slightly deviates from the case of 
spherical symmetry. 
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